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In this work we revisit the definition of Matter Trapping Surfaces (MTS) introduced in previous 
investigations and show how it can be expressed in the so-called dual null formalism developed for 
Trapping Horizons (TH). With the aim of unifying both approaches, we construct a 2-1-2 threading 
from the 1-1-3 flow, and thus isolate one prefered spatial direction, that allows straightforward 
translation into a dual nul subbasis, and to deduce the geometric apparatus that follows. We 
remain as general as possible, reverting to spherical symmetry only when needed, and express the 
MTS conditions in terms of 2-expansion of the flow, then in purely geometric form of the dual null 
expansions. The Raychadhuri equations that describe both MTS and TH are written and interpreted 
using the previously defined gTOV (generalized Tolman-Oppenheimer-Volkov) functional introduced 
in previous work. Further using the Misner-Sharp mass and its previous perfect fluid definition, we 
relate the spatial 2-expansion to the fluid pressure, density and acceleration. The Raychaudhuri 
equations also allows us to define the MTS dynamic condition with first order differentials so the 
MTS conditions are now shown to be all first order differentials. This unified formalism allows one 
to realise that the MTS can only exist in normal regions, and so it can exist only between black hole 
horizons and cosmological horizons. Finally we obtain a relation yielding the sign, on a TH, of the 
non-vanishing null expansion which determines the nature of the TH from fluid content, and flow 
characteristics. The 2-1-2 unified formalism here investigated thus proves a powerful tool to reveal, 
in the future extensions, more of the very rich and subtle relations between MTS and TH. 
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I. INTRODUCTION 

The study of cosmological expansion in local self- 
gravitating structures is an ancient question in General 
Relativity, dating back to the works of McVittie [1] and 
Einstein and Straus [2, 3]. The difference between these 
two pioneering approaches lies in that the first adopts a 
smooth metric with appropriate limits at local and global 
scales, while the second introduces a junction between 
a local solution and a global one. Since then both ap¬ 
proaches have been explored in the literature [4-12] (see 
also [13-19] and references therein for related issues). 
One alternative approach based on smooth solution was 
explored in the works of [20-23], where, in spherical sym¬ 
metry, the concepts of Matter Trapping Shells or Sur¬ 
faces (MTS hereafter), also referred to as “separating 
shells”, was introduced. This concept, first discussed 
there in perfect fluids, appears as a candidate to divide 
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the scales of global physics from those of local physics in 
self-gravitating systems with a cosmological background. 
The introduction in this framework of imperfect fluids, 
with anisotropic stress and heat flux terms, was later 
proposed in [22, 23]. 

This approach relies on the idea that one should search 
within general relativistic gravity for a transition between 
the cosmological realm characterized by the expansion, 
and, say, an astrophysical domain, where manifestly the 
non-linear collapse of overdensities either produces viri- 
alized structures or more extreme collapsed objects. This 
goal is parallel, and somewhat complementary to under¬ 
standing why, on smaller scales, structures and physics 
altogether seem to be immune to the overall expansion 
of the universe. However, it does not necessarily require 
that the physics on the smaller scales be newtonian, al¬ 
beit it does not exclude it. Black holes and other extreme 
objects may very well be the outcome of the divide that 
we envisage. And we should also refer to other inter¬ 
esting possibilities such as the cracking phenomenon put 
forward by Herrera and collaborators, which is associated 
to anisotropic instabilities of compact objects in Astro¬ 
physics [24-27] 

Indeed, a different aspect of collapsing spacetimes solu¬ 
tions is the formation of dynamical black holes, for which 
there is a robust formalism first proposed in the seminal 
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paper of Hayward [28] and further progress in the domain 
have been steady since then (see [29-31] and references 
there in). Dynamical black holes are identified and de¬ 
scribed by the behaviour of the trapping horizons (TH). 
It is therefore of great interest to keep this viewpoint in 
mind and explore a formal approach that might share 
concepts. 

In this work we aim at making use of the tools that have 
been proved useful in the study of dynamical black holes 
in order to study the properties of MTS’s in collapsing 
spacetimes with a cosmological background. We write 
the equations that define the MTS and those defining the 
TH in terms of the same scalars and operators, by finding 
the relationship between the scalars in the 1-1-3 formalism 
— which have been used in the former literature in the 
MTS — and the scalars in the 2-1-2 formalism — which is 
mainly used in the TH literature. We then present both 
points of view in the common 2-1-2 formalism. We are 
led to obtain a form of the MTS conditions which avoids 
derivatives of the 2 -expansions, and that thus establishes 
algebraic equations relating the 2 -expansions, the curva¬ 
ture of the 2 -spheres ( the term) and the energy mo¬ 
mentum tensor. This illustrates how the consideration of 
the role of matter within the present geometrodynamical 
formalism may be of importance. 

Finally we outline further progresses in interpretations 
and then conclude. 


II. THEORETICAL FRAMEWORK 


We first recall the main concepts and conditions defin¬ 
ing both MTS and TH, respectively starting from a 1-1-3 
framework and from a 2 - 1-2 framework and proposing 
some clues about translation of these in the respective 
framework of each other’s concept. 


A. Definition of MTS in 1+3 formalism 


1. Local definition 


In [20-23], a formalism based on a mixture of the 3+1 
(ADM [32]) metric decomposition and the 1+3 [33-36] 
spacetime decomposition was adopted with the purpose 
of separating local from global physics. This formalism 
allows one to split quantities along and orthogonal to the 
direction of the fluid flow. Following Lasky and Lun [37- 
40] , the metric follows the form referred to as generalised 
Painleve-Gullstrand (GPG) [37, 41] 

Here the functions of the temporal and radial coordinates 
t and R that we denote by a,l3,E and r, respectively, 


are the lapse and radial component of the shift vector 
(from the nomenclature of [32]), the energy-curvature of 
the frame-flow-orthogonal hypersurfaces, and the areal 
radius. Also H represents the solid angle. For a per¬ 
fect fluid, the MTS are defined by two equations [20], 
establishing distinct conditions. The first, a kinematic 
condition, guarantees the instantaneous conservation of 
the spherically symmetric Misner-Sharp mass [42] inte¬ 
rior to the shell in terms of gauge invariant quantities, 
namely of the flow’s optical scalars. 
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where m = is the source fluid flow 4-vector. We have 
introduced the new notation ©3 for its expansion, a de¬ 
notes its shear scalar as defined in [37, 41], £„ denotes 
the Lie derivative along u and r is the areal radius of the 
MTS. Since = u°'dar, and the 1-form Ua = dar is the 
gradient of r, Eq. (2) can be interpreted as the locus A (a 
hypertube) where the flow is tangent to a hypersurface 
foliated by spheres of constant areal radius {u°^na = 0 ). 
In other words, the hypertube defined by Eq. (2) is fo¬ 
liated by the spheres where the area of the fluid shell 
reaches an extremal value. 

The second condition is dynamical, and related to a 
generalisation of the hydrostatic equilibrium equation 
of Tolman-Oppenheimer-Volkoff (TOV) through a func¬ 
tional that was termed gTOV [20]. It reduces to 


gTOV _ 





in terms of gauge invariants from the flow (or the locus 2 ) 
where u°'da {u^nt) = u^n^a = 0). Eq. (3) reduces to the 
TOV condition for perfect fluids when Eq. (2) is satisfied 
on the same shell. In contrast with the use of the TOV 
equation in stellar-type solutions, where it holds for all 
shells with matter, in the present context it only holds on 
individual shells where Eqs. (3) and (2) are satisfied. The 
former acts as an equilibrium condition, and is completed 
with Eq. (2) as a kinematic constraint to keep the shells 
from escaping equilibrium. If the two conditions, (2) and 
(3) are satisfied simultaneously, we have a dividing shell 
that separates locally an inner region that can collapse 
from an outer region that can expand. Adding initial and 
boundary conditions of cosmological expansion, it is ar¬ 
gued that such division is global. Such dividing shell was 
termed MTS and examples of spherically symmetric so¬ 
lutions can be found in the references [20-23]. One such 
example is the simplest spherically symmetric model of 
dust with a cosmological constant, using a cosmologically 
inspired initial density profile and a Hubble flow to de¬ 
termine the spacetime completely. As an illustration, we 
present some flow lines for such model in Eig. 1, illus¬ 
trating the hypertubes A and 2 ) as well as the flow line 
defining the MTS in that case. 
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hand, as in the examples of those previous papers, con¬ 
ditions b fulfilled at time infinity define the kind of sur¬ 
face, traced back with their integrable properties, with 
the properties we wish to promote as MTS. For the pur¬ 
pose of this paper we therefore propose to use the follow¬ 
ing definition: 

c. Asymptotic Global Matter Trapping Sur¬ 
face Definition: A Matter Trapping Surface is 
defined by a surface that evolves along flow lines 
and verifies condition b at time infinity. 


B. Trapping Horizons in Dual Null formalism 


Figure 1. A sketch of the worldlines for the second dust-|-A 
example of Ref. [20]. The initial condition are typical of cos¬ 
mology. In this case, the points of ^ are the areal radius 
turnaround events while those of T> are its zero acceleration 
events. Both asymptote to the MTS. 

2. Global definition 

In previous works [20-23], the local definition was used 
to characterise surfaces that would globally separate ex¬ 
pansion from collapse with the crucial assumptions that 
some FLRW model would match at spatial infinity and 
that initial velocities should all be directed outwards in 
a Hubble-like fashion, as expressed by the requirement of 
cosmological initial and boundary conditions. 

a. Cosmological boundary conditions: A model 
is declared to have cosmological boundary condi¬ 
tions if its initial velocities are all directed outwards 
and it tends to a given FLRW model at its spatial 
infinity. 

b. Local Matter Trapping Surface Definition: A 

2-surface is locally called Matter Trapping Surface 
iff it verifies Eqs. (2) and (3) simultaneously. 

What set of events would, with those conditions b, consti¬ 
tute a global definition of an MTS was left open. Aesthet- 
ical considerations would drive one towards a continuous 
3-volume made of 2-surfaces over a nonzero mesurable 
set of flow trajectories, so that the local conditions b 
should happen for continuous extended time. However, 
in examples presented in these papers, such requirement 
is not fulfilled. Instead, as seen in the A and dust models 
of [20, 21], and in the radiation with dust model of [22], 
conditions b were fulfilled only for a 2-surface of events at 
time infinity, and the surface defined and identified only 
thanks to integrable properties of these models. 

In retrospect, this is not too surprising: the two condi¬ 
tions b describe necessary properties of the MTS in the 
shape of 3-world sheets. However they are not expected 
to coincide in a nonzero mesurable set. On the other 


A very useful notion for the study of dynamical space- 
times and characterization of black holes in such setups 
was proposed by Hayward in [28], where the Trapping 
Horizons(TH) where first defined. 

The main idea is that, if the spacetime can be foliated 
in codimension-two spacelike surfaces Es,t where s and 
t are the parameters that identify each surface, then we 
can describe the Lorentzian two-dimensional submanifold 
orthogonal to the E’s in terms of dual null coordinates. 
Therefore, each E surface can be understood as a congru¬ 
ence of null curves, in two independent null directions. 
We can then classify the 2-surfaces according with the ex¬ 
pansion of the null curves across the surface, which we’ll 
call here the 2-expansions, to distinguish from the ex¬ 
pansion ©3 of the flow vector across its orthogonal three 
dimensional space. The 2-expansions of the congruences 
defined by a given surface E and a vector field are 
given by 

©(„) = , (4) 

where the 2-expansion of the congruence v is noted ©(.;,] 
and is defined using the projector onto the surface E or¬ 
thogonal to the dual null directions, say and gen¬ 
erally set up as 

^ab = 9ab T ^(a^b) ■ (b) 

In the case of spherical symmetry, we can use the areal 
radius scalar, that is constant on each sphere, to write the 
two-expansions of a vector field v‘^ across the spheres: 

= ( 6 ) 

Choosing two linearly independent null vector fields in 
the subspace orthogonal to the 2-spheres of symmetry, 
and we can define three types of spheres according 
to the value of the null two-expansions across it: 

• If 0(fc)0(;) < 0 on the sphere, it is called a normal 
or untrapped surface. If a region of the spacetime 
can be foliated in normal surfaces, we call it a nor¬ 
mal region, 
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• If 0(fe)0(;) > 0 on the sphere, it is called a trapped 
surface. Depending on the sign of the expansions, 
a trapped surface can be: 

— future trapped or simply trapped if 

0 (fc) , 0(0 < 0 ; 

— past trapped or anti-trapped if 0(fc) , 0(j) > 0. 

If a region can be foliated in terms of trapped sur¬ 
faces we call it a trapped region. Trapped regions 
can also be classified as future or past. 

• If 0(fe)0(/) = 0, then it is called a marginal surface. 
Without loss of generalily, we assume that 0(fc) =0 
on a certain marginal surface. Thus, this surface 
can still be either 


choice of basis in the Lorentzian 2-space orthogonal to 
the spheres. We can chose a pair timelike/spacelike vec¬ 
tor as a basis, or two null linearly independent vectors 
as our basis, one is easily derived from the other. The 
passage from 1-1-3 to a codimension-two formalism is the 
subtlest step in our work. 

We will first introduce the foliation most natural to spher¬ 
ical symmetry, then the 1-1-3 and 2-1-2 projectors that 
allow to separate the flow from its orthogonal hypersur¬ 
faces as well as the symmetry sheet. We will then recall 
the construction of the flow optical scalars before using 
them in an other look at the MTS definition. 


1. Foliations in spherical symmetry 


— outer, if CiQ(^k) < 0; 

— inner, if LiQ^k) > 0; 

— degenerated if = 0; 

A trapping horizon is defined as a Lorentzian 3- 
hypersurface foliated by marginal surfaces. Trapping 
horizons can be boundaries between regions of different 
kind. We can classify trapping horizons as past/future 
and outer/inner in the same way as marginal surfaces. 

In static solutions, TH’s coincides with event horizons 
(EH). But in dynamical solutions, TH and EV corre¬ 
spond to different hypersurfaces in general. Since the 
EV definition is global and the determination of the EV 
is only possible after integration of the whole future of 
the spacetime, in the dynamical case the EV has been 
replaced by the TH as a definition or characterization 
of black holes. Black Holes, in the sense of a very com¬ 
pact object that is created as the result of gravitational 
colapse, are associated with future outer trapping hori¬ 
zons. On the other hand, cosmological horizons are past 
inner trapping horizons in this formalism. For further 
detail on Black Holes dynamics and TH theory we refer 
the reader to [13, 28-31, 43-45] 


III. TRANSLATING 1+3 TO 2+2 FORMALISMS 

A. MTS definition adapted to a codimension two 
foliation 


In the following discussion, we will expand the framework 
of MTS presented in Sec. H A and express its definition in 
a codimension-two formalism, i.e., in a formalism where 
consider the evolution of scalars defined on the leaves of 
codimension-two foliation, instead of the usual formal¬ 
ism, which defines its scalars in a 3-space orthogonal to 
the flow, being a codimension-one formalism. 

We refer to the codimension-one formalism as the 1+3 
formalism, while the codimension-two formalism is re¬ 
ferred to as the 2+2 formalism, but depends on the 


Since, in spherical symmetry, spacetime is naturally fo¬ 
liated in 2-spheres, that are maximally symmetric, it 
is fruitful to adopt a formalism which already explicit 
such symmetry, that is where the tangent space to the 
2-spheres is separated from their orthogonal space. Fol¬ 
lowing Clarkson [46-48], we consider an orthonormal 
basis {■«“, e“, s/ 2 )}, where s/^) and s/ 2 ) form any 

basis tangent to the 2-spheres, is the fluid flow 4- 
vector and e“ is a normalised spacelike vector such that 
M“ea = s“ea = 0, where s“ is any vector tangent to the 
symmetry 2-spheres. However, this decomposition can 
apply to any spacetime for which the flow-orthogonal 
hypersurfaces can be decomposed between a privileged 
vector field and its codimension-two orthogonal hyper¬ 
surfaces. 

Also, considering its equal suitability to the symmetry 
of the problem and usefulness in the study of spherically 
symmetric solutions as well as in the formalism of trap¬ 
ping horizons, we shall use the 2+2 formalism, by using 
the basis vector from the above to build a null basis of two 
future directed and linearly independent vectors, {fc“, 1“}. 
Again, the codimension-two hypersurfaces orthogonal to 
that basis need not, in general, be spherically symmetric. 

We restrict ourselves to the case where the sheets are al¬ 
ways spacelike, that is, all their tangent vectors are space¬ 
like. Here we define the notation and relations between 
basis vectors and projectors that shall be of foremost im¬ 
portance for this work. 

Starting from the normalised and orthogonal vector fields 
{u°‘,e°‘}, satisfying: 


and e“ as 


= -1 , 

(7) 

e“ea = 1, 

(8) 

u°-ea = 0. 

(9) 

nding dual null basis {A:“ 

, /“ } using 

/c“ = + e“ , 

(10) 

- e“ , 

(11) 
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that satisfies, by construction 

fc“fca = na=0, (12) 

k^la = -2. (13) 

Imposing the constraints (7), (8) and (9) does not defines 
our basis uniquely, since we can build a new pair of basis 
vectors by a local Lorentz boost 

= cosh O' + sinhw e“ 
e'“ = cosh a; e“ + sinh w , (14) 

where w = w(a;“), and get a new orthonormal basis 
{u'“, e'“} satisfying the same constraints. In the corre¬ 
sponding null basis, the boost translates into 

, 

(15) 

We have then defined two sets of bases for the orthog¬ 
onal space of codimension-two foliation, one set being 
orthonormal, the other is dual null. We have explicited 
the constraints and freedoms of redefinition such basis 
should enjoy. In this work, we shall choose u°‘ in the di¬ 
rection of the matter flow, which also sets e“ from the 
foliation up to a sign and, subsequently, the basis null 
vectors without ambiguity. 


2. Projectors 


In this section, although we shall remain in the 4- 
dimensional case, we will strive to keep the notation in¬ 
dependent of dimensionality as much as possible. 

As is customary, for instance in 1-1-3 [33-36], the projec¬ 
tor onto the subspace orthogonal to reads 

hah — Qah 4 “ UaUh 7 (^^) 

and, analogously as in [46-48] or in [28, 29, 49], the fur¬ 
ther projector onto the 2-subspace orthogonal to u°‘ and 
e“ or in and 1°“ respectively writes 


A^o6 — 3ah 4“ UaUb ^a^h — 3o,h 4- k^a^b) • (^^) 

In this framework, we define the notations for contrac¬ 
tions of general rank-2 tensors Aab with respect to these 
projectors: 


Ah = h‘^^Aab; AN = N‘^^Aab. (18) 

We are interested in projecting physically meaningful 
tensorial quantities (in particular the covariant derivative 
of the flow) between their tangential and orthogonal com¬ 
ponents to the flow. However, we are also interested in 
decomposing the orthogonal components into their trace 


and traceless parts. ^ . This is achieved with the def¬ 
inition of the projector on the subspace of symmetric 
trace free rank 2 tensors on the 3-space, which includes 
the tracefree operator that leaves unaffected the leaves 
of the 2-foliation normal to the dual null subspace. We 
name this latter operator the STFPr {symmetric trace 
free projector):^ 


Pab — hab h^CaCb — P^c^a^b ■ (^^) 

Following its definition, we obtain the symmetric trace- 
free 3-spatial component along that STFPr of a rank 2 
tensor from the contraction: 

pab A 

(20) 

where, in 4-dimensions, we get 

p2 ^ pabp^^ ^ ( 21 ) 


Finally, we can show an important identity for the trans¬ 
lation between scalars in 1-1-3 to 2-1-2 formalism. We first 
note Ae = Aabe-°‘ef’ the doubly projected component of 
Aab along e“. Applying the STFPr to an arbitrary ten¬ 
sor Aob, thanks to Eqs. (19) and (20) we get the following 
identities 


P’^^Aab = P^a = Ah- KAa = An- NfAe . ( 22 ) 

Eliminating Ag and solving for A^v, we obtain 

p2 

An = + -j^a , (23) 

and, using (21), we finally get 


An _ ^ 
~ 


(24) 


3. Flow scalars 


With the aim to study MTS’s, we need to investigate 
the geometric properties of the flow as described by 
the independent components of its covariant derivative 


^ In the particular case of spherical symmetry, they can both define 
scalars and the only degree of freedom orthogonal to the flow is 
along the radial e“, the spacelike normal to the 2-spheres. 

^ In [46—48], it can be seen as the traceless combination of the 2- 
projector with the projector along the spacelike normal to the 
leaves of the foliation. The interest of this new projector lies in 
the case of spherical symmetry, in which we shall restrict, where 
all trace free symmetric tensors can be written as a{x^)Pab, 
where a is a scalar function. In other words this representation 
is one dimensional. 
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VaUh- Our goal is to define physically meaningful scalars 
adapted to our formalism. 

Specializing in the 4-dimensional case, in the 1-1-3 formal¬ 
ism, the flow is decomposed as 


VoUb 


, 03 , 
-UbUa H- 


(25) 


derivatives is proportional to the induced metric on the 
2 -spheres Nab, a consequence of the maximal symmetry 
of the spheres. In other words the only non-vanishing 
component of the flow derivatives is the trace component 
(for things living on the 2 -spheres) and thus its traceless 
parts are vanishing. 

Using the definitions of Sec. Ill A 2, and setting Aab = 
VaUb, we realise 


where ©3 = h‘^^VaUb can be recognised as the volume ex¬ 
pansion scalar of the fluid, that we shall also call as the 
3-expansion, to distinguish it from the 2-expansion, de¬ 
fined below; aab is the shear tensor and uJab is its vorticity, 
corresponding to its projected antisymmetric component. 
In spherical symmetry, there is no vorticity in finite- 
streamed flows, thus all symmetric trace-free tensors 
living on the 3-space orthogonal to the flow are pro¬ 
portional to Pab, which implies aab = (^{x°^)Pab and 
we shall redefine the shear scalar as this a{x°‘) from 

cr(x®) = = ^P°‘^'VaUb (recall the classic defi¬ 

nition [33-36] noted here aE = \cr°-^aab = V^cr). 

However, comparing with the results of Sec. Ill A 2 with 
the specialisation Aab = NaUb, we obtain 


Ah = Q 3 , a = a (26) 


An = ©(«)■ 


(30) 


Hence, using the identity (24) and Eqs. (26) and (30), we 
obtain 



(31) 


The above result is the fundamental relation we need to 
bring MTS’s and TH’s to a unified formalism. 


Defining Matter Trapping Surfaces in 2+2 

In the 1-1-3 language, MTSs are characterized by the sys¬ 
tem of Eqs. (2) and (3). Erom Eq. (31), we can express 
them in terms of 2 -expansions as 


The next step is to write the optical scalars of the flow 
in the codimension-two formalism. Erom now on we will 
explicitly suppose spherical symmetry. 

The orthogonal subspace is generated by the basis 
{m“, e“}. We can define 2-expansions, describing the rel¬ 
ative variation to the sphere’s area ( 2 -volume or surface) 
following an infinitesimal displacement along each con¬ 
gruence orthogonal to the sphere (i.e. u°‘ and e“). As 
the spheres have codimension-two, the 2 -expansions of 
the basis vector are of natural interest: 

©(„) = , 

0(e) = N^^Vaeb ■ (27) 

We can then construct the mean curvature vector as 

/C“ = -©(„)«“+ 0 (e)e“ = ^V“r, (28) 

such as, for any in the orthogonal subspace {u, e), we 
have the 2 -expansion along defined as 0 („) = v^K-b- 
In particular, for the 2-|-2 basis vectors, we find 


0(.)=O, (32) 

^-00) + ^=O. (33) 

This is a surprisingly simple form for the conditions, that 
suggests that the 2 - 1-2 formalism probably best captures 
the symmetries of the problem, only requiring derivatives 
in the flow direction m“. Furthermore, we wish to express 
these equations in a dual null 2 - 1-2 basis such that we can 
use the dual null formalism usual in black hole dynamics 
and TH’s for MTS’s. Equations (32) and (33) in terms 
of /c“ and P expansions yield 

0(fe) + 0p) = 0 , (34) 

/:fc0(fc) + + CiO^k) + ^ ^ = 0 , 

(35) 

the equations of MTS in the 2-1-2 formalism. Moreover, 
we point out that the left hand side of Eq. (35) is equiv¬ 
alent to —8 sdUX, 

r 

B. Relation to Matter Sources 


0(fe) k Ka 0(u) + 0(e) , (pjjg 2-1-2 projection of EFE can be decomposed into 

0(i) = PlCa = 0(u) — 0(e) • (29) scalar equations along the null vector 


In spherical symmetry there is neither 2-shear nor 2- 
vorticity as the only non-vanishing component of the flow 


^fc 0 (fc) — i^(fe) 0 (fe) + 



-SnTabk^^k^ , 


(36) 
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with the normalised cross projection of its covariant 
derivative 

m , (37) 

along the null vector 

02 

A0(i) - + -^= -8TTTabn ^, (38) 

with the normalised cross projection of its covariant 
derivative, 

^^(0 = (39) 

the latter two equations being usually referred to as the 
null Raychaudhuri equations, and the mixed projection 
of the symmetric EFE that yields 

2 {^kQ{l) + -C/0(fc) + 0(;)7i + 0(fc)7(fc)) + 0(fc)0(Z) 

-^^^^=87rTabl^k\ (40) 

where the linear coefficients of the expansions are related 
to the dual null commutators, their Lie derivatives (recall 
[kj]‘^ = Ckn, 

7(fe) = ^?a[A:,;r, (41) 

7(0 = (42) 

We can relate the coefficients (37), (39), (41) and (42) in 
useful a way. Let’s consider the difference — 7(fe): 

I'm - 1(1=) = ^ [k-fV.k, - l.k'-Vtl- + l.fVtk"] = 

-j^l(k"l- + n-) V.t,] = ^ = 

nj t Q hjQ L 

-V(i) + 7(Z) , (43) 

where we used the fact that Va{k^lb) = 0 by construction, 
and the fact that the expressions are all symmetric in 
and 1°“. Writing Eq. (43) in terms of the space and time 
orthonormal vector basis, we obtain: 

2 

V{k) - l{k) = {u^'^aSb - e^’VaUb) = 

— (-2e''V,n,) = -j^e^Ub = 2A, (44) 

rh Lq r\j L (2 

where = u'^VaU^ and we have used the fact that 
VaU^^b = 0, by construction, and the fact that we have, 
in spherical symmetry, 

(e°Uc) = Tle^, (45) 

so |7l| = \/ U°‘Ua ■ (46) 


Summing all Raychaudhuri equations, Eqs. (36), (38) and 
twice (40), one can isolate the Dynamical condition (35) 
and the Kinematic condition (34) of the MTS showing the 
MTS to remain with no dynamics, as expressed in the TH 
2+2 formalism. The summed Raychaudhuri equations 
can be written thus 


>Cfe0(fe) + + CkQ(i) + CiQ(^k) 


(0(fc) + 0(/))' 


+ 0(0 {li - + 0(fc) (7(fc) - ’^(k)) + 0(fe)0(O “ 2 


k°-la 


= -SirTab 




<+■ 


•^fc0(fe) + + CkQ{i) + CiQ(^k) + 


(0(fe) + 0(;)) 
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+ — + 87rTaf,4e“e^ — —0(fe)0(i) + 271 (0(fe) — 0(o) 

= 0(fc) + 47l0(fe) - [27l + 0(fc)] (0(fc) + 0(i)) , (47) 


the latter equation using the definitions (10), (11), (7), 
(8) and (9) in the curvature and source terms. Note that 
the source term corresponds, in spherical symmetry, to 
the isotropic pressure P = Tabe°'e^, while the spherical 
curvature has a definite sign. The left hand sign term 
between square brackets presents the vanishing part of 
the dynamics condition (35) S, while the last parenthesis 
presents the vanishing part of the kinematic condition 
(34) .S. Note also that the symmetry between k and I 
has to be broken in the writing of Eq. (47) where we 
isolate Q(k) but that this symmetry remains underlying. 
This combination of the EFE can be used to propose an 
alternative gauge invariant 2+2 definition for the gTOV 
functional in (3) as 


gTOV_ 1 0(fc)0(o 


2r‘ 


- AnTabe^e^ - 

A 


7 (0(fc) - 0(i)) ■ (48) 


At this point it is worthwhile to remind the definition of 
the Misner-Sharp mass M as [42, 43] 

g^^^dardbr = 1 - — , (49) 

r 

which, considering Eq. (28) and Eq. (29) , allow us to 
relate the Misner-Sharp mass with the null 2-expansions: 

/C“/Ca = -0fc0, = 4 fl-—) • (50) 

\ r J 


We can now point out that the reader can recognise the 
perfect fluid gTOV of [20], 


gTOV = 


—— P' + 47rPr + 
p + P 


M' 


(51) 


















where p is the energy density of the fluid defined as p = 
T°'^UaUb and the' stands for radial derivation, in Eq. (48). 
The right hand side of Eq. (48) displays the pressure term 
in its second term, the mass term in its first and third 
terms leaving the pressure gradient in Eq. (51) for the 
last term, that can be identified as 


^ ^ ^ u' 

-^0(e) — -- X ^ _j_ pP ■ 


(52) 


IV. DISCUSSION AND CONCLUSIONS 


The form of the MTS conditions given in Eqs. (47) 
and (53) has the advantage that the derivatives of the 
2-expansions are eliminated, and what remains are alge¬ 
braic equations relating the 2-expansions, the curvature 
of the 2-spheres ( the term) and the energy momen¬ 
tum tensor. In terms of coordinates, we are left with 
a system of first order nonlinear differential equations, 
which means that we greatly simplify the computation of 
the .ft, 2) and MTS surfaces in a given solution or initial 
condition. 

If we restrict ourselves to the local MTS definition, con¬ 
dition b of Sec. II , where (35) and (34) apply together, 
the dynamical evolution Eq. (47) becomes (recall then 
0(0 = - 0 ( fe )) 

+ 4M0(fc) =4 + SirTab . (53) 

In this case the MTS equation reduces to just one first 
order equation, as a kinematic condition. 

In addition, the 2-|-2 formalism renders explicit that the 
condition (34) can only be fulfilled in the so-called normal 
regions of spacetime, which are regions foliated in normal 
surfaces, or on extremal horizons, where 

0(fe)=0(q=O. (54) 

This property makes sense intuitively since trapped re¬ 
gions are related to collapse (future trapped) and to cos¬ 
mological expansion (past trapped). Furthermore any 
MTS cannot belong to either regions. Extremal horizons 
can be understood as situations when two trapping hori¬ 
zons coincide. In this case, the extremal horizon itself 
gives a separation between collapse and expansion. In 
other words, the MTS has to be contained between inte¬ 
rior THs and an exterior TH, i.e. between black holes’ 
and cosmological horizons, as expected. 

On the other hand, we can use the tools we developed 
to relate the TH’s to the MTS quantities, in the same 
framework. In fact for a TH characterised by 0(fe) = 0, 
and using the version (35) of (3), the dynamics simplifies 
into 


gTOV 2 
M0(q =4^-- 


(55) 


in addition to the relevant evolution equations, Eqs. (38), 
in that case. Moreover, Eq. (55) is interesting, because 
it shows that the gTOV functional appears as a balance 
term against the curvature term, and the pressure term 
to give the sign of 0(q, once the sign of the acceleration 
is known. We are provided with a relationship between 
two quantities that convey a notion of contraction or ex¬ 
pansion, the gTOV functional and the null 2-expansion 
0(q at a TH. 

We conclude emphasizing that this formalism suggests 
the possibility to study more subtle relations between 
the trapping horizons and MTSs to describe the relevant 
quantities on an equal footing. It can also serve as a basis 
for an MTS definition independent of spherical symme¬ 
try, in the style of the formalism developed to study dy¬ 
namical black holes proposed in [28]. Our present results, 
which reveal the important role of the characterization of 
matter within the present formalism, find a follow up in 
the further investigation of specific matter contents. The 
latter was beyond the bare scope of the present work, and 
will be addressed in future work to deepen our under¬ 
standing of the dynamical implications of the interplay 
between matter and geometry. 
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Appendix A: Derivation of covariant focusing 
equations 


The form of the Raychaudhuri equation (36) and (38) are 
obtained from the usual form found, for example, in [50] , 

02 

•^if0(if) =-^ ~ <^{K)ab<^(K) + ^{K)ab^(K) ~ RabK°' 

(Al) 

by imposing spherical symmetry, which eliminates the 
2-shear and 2-vorticity terms, and relaxing the affinity 
condition on AT^, K°‘VaK^ = 0. Gonsidering a colinear 
but non-affine null vector field and replacing 

it by AT® in Eq. (Al), we obtain 

02 

21fe0(fc) ~ H-+ Rabk°'k^ = 0 , (A2) 






with = k^^daV. Replacing K'^ by fc“ in the affinity 
condition K°'VaK^ = 0, we obtain 

Vafc'’ = k'^daVk^ = v(k) ■ (A3) 

We can use the same reasoning on the /“ Raychaudhuri 
equation in order to obtain Eqs. (37) and (39). 

The mixed kl projection of EFE can also be obtained 
starting by its expression in a dual null coordinate form 
obtained, for example, in Ref. [51]: 

2 ('^-®(+) + ^+®(-)) = “®(+)®(-) + + 87rT+_ , 

(A4) 

where the spherically symmetric line element takes the 
form 

ds^ = 2(7-1_du_|_dM_ + r'^{u+, u_)dfl^ . (A5) 

This equation is valid for a null set of coordinates 
{m-i-jU-}, that satisfies [i9-|-,9_] = 0. A general dual 
null basis {fc“, /“} does not commute and does not corre¬ 
spond to the derivatives along any coordinate. Using the 
same reasoning as earlier, we suppose that the basis are 
related in the following way: 

= e''(a_H)“, (A6a) 

P = e’^(a_)“. (A6b) 
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Replacing Eqs. (A6a) and (A6b) in Eq. (A4), we obtain 

i (£fe0p) + CiQ^k)) = -e(fc)0p) + + 87rr„f,rfc''+ 

^{Q(i)k‘'dctp+ Q(k)l''dcv) ■ (A7) 

The last step is to characterize in a covariant manner the 
functions k‘^dc'il> e q'^dcT] in Eq. (A7). Using the commu¬ 
tators: 

0= [a+,5_] = = 

= e-^-^ {[k, l]^ - {k-da^P)l^ + {l^dar])k^} . (A8) 

Finally, we obtain 



(A9) 

ra.i = prl''T^ 

(AlO) 
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